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We search for direct pair production of suporsymmetric top quarks and supersymmetric bottom 
quarks in proton-antiproton collisions at ^ = 1.96 TeV, using 295 pb~^ of data recorded by the 
Collider Detector at Fermilab (CDF 11) experiment. The supersymmetric top (supersymmetric 
bottom) quarks arc selected by reconstructing their decay into a charm (bottom) quark and a 
neutralino, which is assumed to be the lightest supersymmetric particle. The signature of such 
processes is two energetic heavy-flavor jets and missing transverse energy. The number of events that 
pass our selection for each search process is consistent with the expected standard model background. 
By comparing our results to the theoretical production cross sections of the supersymmetric top and 
supersymmetric bottom quarks in the minimal supersymmetric standard model, wo exclude, at a 
95% confidence level in the frame of that model, a supersymmetric top quark mass up to 132 GeV/c^ 
for a neutralino mass of 48 GcV /c? , and a supersymmetric bottom quark mass up to 193 GeV /(? 
for a neutralino mgiss of 40 GeV/c^. 

PACS numbers: 12.60.Jv, 13.85.Rm, 14.80.Ly 
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I. INTRODUCTION 

Supersymmetry (SUSY) is an extension of the stan- 
dard model (SM) of particle physics that overcomes some 
of the theoretical problems in the SM by introducing a 
new degree of freedom [ij. In this model a bosonic su- 
persymmetric partner is assigned to every SM fermion 
helicity state, and a fermionic superpartner to every SM 
boson. Thus, the SM quark helicity states and 
acquire scalar partners qL and q^. The mass eigenstates 
of each supersymmetric quark (squark) can be a mixture 
of their weak eigenstates, quantified by a mixing angle 
9. The difference in the mass eigenvalues depends on 
several factors. In the case of the supersymmetric top 
quark (stop), due to the large top quark mass and the 
large value of its Yukawa coupling constant (Higgs-to- 
top coupling), there can be a significant difference in the 
mass between the two mass eigenstates ti and t2- In the 
case of the supersymmetric bottom quark (sbottom), a 
large mass difference between the two mass eigenstates 
hi and 62 can occur if the ratio of the vacuum expecta- 
tion values of the two Higgs fields expected in SUSY is 
large [2]- In both cases it is likely that the less massive 
stop {ti) and sbottom (61) could be hghter than the first 
two generations of supersymmetric quarks. In fact the 
ti could be lighter than the top quark if the top Yukawa 
coupling strength or the stop mixing is strong enough. 
The requirement of a light stop is also a feature of many 
baryogenesis models jSj]. 

In this paper, we describe two analyses searching for 
stop and sbottom production at the Tevatron using a 
data sample of 295 pb^^ integrated luminosity collected 
with the CDF II detector. Both analyses are performed 
within the i?-Parity [i] conserving minimal supersym- 
metric standard model framework. A consequence of R- 
Parity conservation is that all SUSY particles are pair 
produced, and the lightest SUSY particle (LSP) is sta- 
ble. If the LSP interacts weakly it is a good candidate 
for cold dark matter and it escapes our detection. 

At the Tevatron, stop and sbottom are expected to be 
produced in pairs mainly via gg fusion and qq annihila- 
tion, as shown in Fig. [T] At leading order, their pro- 
duction cross sections depend essentially only on their 
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masses. For a center-of-mass energy of 1.96 TeV, the 
next-to-leading order (NLO) cross section, calculated 
with PROSPINO [5, ranges from 50.3 pb to 0.25 pb 
for stop and sbottom masses from 80 GeV/c^ to 200 
GeV/c^. In the calculation, the renormalization and fac- 
torization scales are set to the mass of stop or sbottom 
(Qj-f — ) , and the parton distribution functions 

(PDFs) are from CTEQ6m ^] . 




FIG. 1: Leading-order Feynman diagrams for pair production 
of stop and sbottom at the Tevatron. 

The stop and sbottom can decay in many channels, 
depending on the mass difference between stop/sbottom 
and other SUSY and SM particles. Here we consider the 
SUSY parameter space where the stop and sbottom are 
relatively light. In the case of stop, the flavor changing 
loop decay ti —> cxi dominates if m^^ > nic + m-^o , but 
TOjj^ < mf) + Tn^± and m^^ < m\\r + rrif, -|- m^^o. In the 

sbottom search, the 61 — > 6x1 is the only relevant decay 
if m^^ > TOfc -I- m-j^o , but < rrib + m^o and mg^ < 
TT^t + "^j^± ■ The neutralinos {xi 234) ^^'^ the charginos 
{X12) the SUSY partners of the electroweak bosons 
and are labelled in order of increasing mass. Therefore, 
we search for the processes pp titi (cx?)(cxi), 
and pp bibi (^Xi)(^Xi); shown in Fig. (5] We 
assume Xi is the LSP. Thus, the experimental signature 
for stop and sbottom pair production processes is a pair 
of acoUinear heavy-flavor jets (ie. b and c jets), and large 
missing transverse energy {I^t ) coming from the escaping 
LSPs. The searches assume that both stop and sbottom 
decay very close to the interaction point. For =100 
GeV/c2 (TOg^ = 160 GeV/c^) and m^o = 60 GeV/c^, the 
expected lifetime of stop (sbottom) with only the decay 
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h cx? (&i ^ &Xi) is of the order of 10'^^ (IQ-^S) 
seconds which corresponds to a natural decay length of 
0.3 /xm (3 X 10-9 fim) [3,i]. 



hadronic calorimeters is (t{Et)/Et 



Xi 



W 




FIG. 2: The decay channels of stop and sbottom considered 
in this paper. 

Previous searches for stop and sbottom have been per- 
formed at LEP and at the Tevatron 9, 10, 11, 12, 13, 14, 
[ll,[ll,[l3|. For the search topology studied in the present 
analysis, LEP excludes stop (sbottom) masses smaller 
than « 100 GeV/cP 100 GeV/c^), independent of 
the difference between stop (sbottom) and neutralino Xi 
masses @ . Recent results from the D0 Collaboration in 
the same topology [H, , based on Run II data, have 
extended Tevatron's Run I reach [13, [lH by excluding 
stop masses up to 133 GeV/c^, and sbottom masses 
up to ^ 220 GeV/c^. These results are also shown in 
Fig. 0] and Oof this paper. 



II. THE DETECTOR AND DATA SAMPLE 

CDF II is a general-purpose detector that is described 
in detail elsewhere [l^ • The components relevant to this 
analysis are briefly described here. The charged-particle 
tracking system is closest to the beam pipe, and con- 
sists of multiple layers of silicon micro-strip detectors, 
which cover a pseudorapidity region |r/| < 2, and a large 
open-cell drift chamber covering the pseudorapidity re- 
gion jryl < 1 [iqI. [20|. [2H . The silicon microstrips of the 
silicon detectors have a pitch of 25 to 65 /xm, depend- 
ing on the layer, thus allowing a precise measurement 
of a track's impact parameter with respect to the pri- 
mary vertex. The tracking system is enclosed in a su- 
perconducting solenoid, which in turn is surrounded by 
calorimeters. The calorimeter system 22| is organized 
into electromagnetic and hadronic sections segmented 
in projective tower geometry, and covers the region 
\ri\ < 3.6. The electromagnetic calorimeters use lead- 
scintillator sampling, whereas the hadron calorimeters 
use iron-scintillator sampling construction. The trans- 
verse energy resolution of the electromagnetic calorime- 
ters is a{ET)/ET = I — ®2% for the central region 



(M < 1), and a(i?)/£; = 
region (|?7| > 1). The transverse energy resolution of the 



®3% 

80"''" ® 5% for 



1% for the forward 



for the central region, and a{E)lE = ^=7„ 

the forward region. The present analysis exploits the in- 
formation of the central muon system, which is located 
outside of the calorimeter and covers the range |?7| < 1. 

The data sample for this analysis was collected using a 
Ifx +jets trigger, which is implemented in three levels of 
online event selection. The I^t is defined as the energy 
imbalance in the plane transverse to the beam direction 
dH, and a jet is defined as a locahzed energy deposition 
in the calorimeter. In the first and second levels of the 
trigger, Ifx is required to be greater than 25 GeV and 
is calculated by summing over calorimeter trigger towers 
with transverse energies above 1 GeV. In the second level 
there must be at least two jets with Et > 10 GeV. In the 
third level, I^t is recalculated using the full calorimeter 
segmentation with a tower energy threshold of 100 MeV 
and is required to be greater than 35 GeV. 

In the oSline processing, jets are reconstructed from 
the calorimeter towers using a cone algorithm with fixed 
radius AR = A-rf- + A^^ ^ 9 4 -^^ rj-cj) space [Hi. The 
jet Et measurements and I^x £^re corrected for detector 
effects di]. 

A fraction of events passing the trigger is not from pp 
collisions, but from beam halo and cosmic ray sources. 
To remove these events we examine the event electromag- 
netic fraction Fem and charged fraction Fch ■ The Fem is 
the ratio of the energy measured by the electromagnetic 
calorimeter to the total energy contained in jets of cone 
radius AR = 0.4 with Ft > 10 GeV and \t]\ < 3.6. The 
Fch is the fraction of the jet energy carried by measured 
charged-particle tracks {pt > 0.5 GeV/c) averaged over 
jets with jryl < 0.9. The beam halo travels parallel to the 
beam axis and deposits most of its energy in the elec- 
tromagnetic section of the calorimeter or in the hadronic 
section. By requiring Fem > 0.1 we reject events which 
contain little energy in the electromagnetic section of the 
calorimeter. A cosmic ray traversing the CDF detector 
can deposit energy in the calorimeter without register- 
ing a track in the tracking detectors. If the beam halo 
background events and cosmic ray background events, 
as described above, do not overlap with beam crossing 
events that produce hard pp collisions, then there will 
be little activity in the tracking detectors. Therefore, by 
requiring Fch > 0.1 we reject backgrounds that have lit- 
tle tracking activity. More detailed explanations of these 
two variables and how they reduce the non pp collision 
events are described in [2^. We also reject events if the 
reconstruction cone of any jet in the event enters an unin- 
strumented region of the calorimeter. 



III. BACKGROUND SOURCES AND EVENT 
SELECTION 

The dominant backgrounds to the stop and sbottom 
searches in the jets and I^t signature are production of 
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multi-jet, W oi Z boson with jets, single top, tt, and di- 
boson {WW /W Z / Z Z) final states. We use the ALPGEN 
generator to simulate the W and Z boson plus parton 
production, with herwig used to model parton showers 
[271 [28|. Multi-jet, top quark, and di-boson production 
are simulated with pythia [23]. All the Monte Carlo 
(MC) SM background samples were generated using the 
CTEQSl PDFs [30]. In generating the multi-jet sample, 
we select events where there are 6 or c outgoing par- 
tons/hadrons in the final state [heavy- flavor (HF) multi- 
jet]. To normalize the SM top background samples, we 
use the NLO cross section values for single top and tt pro- 
duction [31I, m, [11] . We use the mcfm program to obtain 
the NLO cross sections for W/Z-\-]Qis and di-boson pro- 
duction [sH]. The HF multi-jet sample is normalized 
using kinematic regions in the data that are dominated 
by multi-jet production. In multi-jet production events, 
the I^T is usually due to jet energy mis-measurement. In 
this case the IfiT tends to point in the same direction as 
the jet whose energy is mis-measured. In the ^T+jets 
data sample used in this analysis, the events at higher 
I^T are dominated by non-multi-jet SM contributions. 
Therefore, the multi-jet dominant kinematic regions are 
low I^T regions (50 < Ifr < 70 GeV), and regions where 
a jet is aligned in the direction of Ifx ■ We obtain an 
average normalization factor kHFmuiu-jet — 1.46 ± 0.37 
where the uncertainty is mainly due to the uncertainty 
from jet energy calibration and resolution [2^. The 
kHFmuiti-jet factor is used to normalize the HF multi-jet 
sample to estimate the HF multi-jet contribution in the 
signal region. 

Data selection is optimized by maximizing the statis- 
tical significance of a simulated stop/sbottom signal over 
the expected background events in the data. The opti- 
mization is performed prior to examining the signal re- 
gions of the data. As the signal production cross sec- 
tion and event kinematics (for example jet Et, Ifr ) 
could vary significantly across stop and sbottom masses, 
we determine separate sets of optimized cuts for three 
mass ranges. The "low", "medium", and "high" mass 
ranges used for the stop (sbottom) search are < 100 
GeV/c2, 100 < < 120 GeV/c^, and m^^ '> 120 

GeV/c2 (m^^ < 140 GeV/c^, 140 < m^^ < 180 GeV/c^, 
and mg^ > 180 GeV/c^). The selection cuts for the stop 
and sbottom are summarized in Table [T| and [TT| respec- 
tively. At the initial stage of event selection, the data 
sample is dominated by multi-jet background events. 
We employ several selection cuts to minimize their con- 
tribution. To reduce multi-jet background and avoid 
regions where the I^t trigger is inefficient, we require 
> 50 GeV. The trigger is - 75% efficient in this 
region. Next we require that there be only two or three 
reconstructed jets in ]77] < 2. Events with any addi- 
tional jets with Et > 8 GeV and 2 < ]77] < 3.6 are 
rejected. Most of the time the two highest-i?T jets in 
multi-jet events are anti-parallel. However, this is not 
the case in the stop (sbottom) pair production since the 
c (6) jet recoils against the Xi in the stop (sbottom) de- 



cay. Therefore, we require that the opening angle be- 
tween the two highest Et jets be less than 160 degrees 
in the plane perpendicular to the beam. The large I^t in 
multi-jet events that survive the earlier cuts is usually 
due to jet energy mis- measurement. Thus, to further re- 
duce this background we require a minimum azimuthal 
separation between the direction of the jets and I^t of 
A(f>{Jet, I^t) > 45°. For multi-jet events that pass the 
minimum azimuthal separation requirement, but have 
large I^t due to jet energy mis-measurement, the magni- 
tudes of the I^T and the second leading jet's transverse 
energy are often anti-correlated. Therefore, we require 
the sum i?T(Jet2)-|-^T to be above the values listed in 
Table [J and [ni 

In stop (sbottom) pair production, where each stop 
(sbottom) decays into a charm (bottom) quark and a 
neutralino, the vector sum of the transverse energy of 
the two c jets (5 jets) should balance against the miss- 
ing transverse energy from the escaping neutralinos. We 
explore the correlation between the missing transverse 
energy and the transverse energy of the first and second 
leading jets through the variable E^ , which is de- 

G J O J J12MET ' 

fined as 



E^ 



{Etji+Etj2+1?t)- 



Ej 



Et 



Etji + Etj2 



(1) 



Etji and Et,2 are, respectively, the vectors of the trans- 
verse energy of the first and second leading jets. For pair 
production of stop (sbottom), E^^^^^^^^ is expected to 
be insignificant, i.e. within the calorimeter jet energy 
resolution. For multi-jet production, due to its high jet 
multiplicity, the missing transverse energy caused by the 
mis-measurement of a single jet's energy, may not bal- 
ance against the transverse energy of the first and sec- 
ond leading jets. In top-quark pair production, the miss- 
ing transverse energy caused by the escaping neutrino 
in the W boson decay also does not necessarily balance 
against the first and second leading jets in the event. 
Therefore, E^ is large in these background events. 

To reduce the contribution from top quark production, 
and to further reduce multi-jet background, wc require 
i?T < 15 GeV for both stop and sbottom analyses. 

J-J12MET ^ ^ 

To reduce the background contribution from 
M^/Z-|-jets and top quark production, we reject 
events with one or more identified leptons. Candidate 
electrons must have a track associated with a cluster 
in the electromagnetic calorimeter with Et > 10 GeV. 
Its electromagnetic to hadronic energy ratio and shower 
profile, must also be consistent with that expected for 
electrons. Candidate muons are identified as tracks with 
Pt > 10 GeV/c that extrapolate to hits in the muon 
chambers and to energy deposited in the calorimeters 
consistent with a minimum ionizing particle. To increase 
lepton detection efficiency, candidate leptons (electrons, 
muons, taus) are also identified by isolated tracks with 
Pt > 10 GeV/c. The number of tracks associated 
with the first and second leading jets should be four or 
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TABLE I: The event selection cuts for pair production of stop in the "low", "medium", and "high" stop mass regions. Jetl, 
Jet2, and Jet3 are, respectively, the first, second, and third leading jets. The cuts that are listed only under the "medium" 
column are common for all three mass ranges. 



Mass range 


Low 


Medium 


High 


(GeV/c^) 


< 100 


100 - 120 


> 120 


(GeV) 


> 50 


> 50 


> 50 


Et Jetl (GeV) 


> 35 


> 45 


> 55 


Et Jet2 (GeV) 


> 15 


> 15 


> 25 


Et Jets (GeV) 


> 15 


> 15 


> 15 


\t) \ : Jetl, Jet2, Jet3 




1 »?i |< 1.2, 1 r)2 \< 1.5, 1 7?3 |< 2.0 




Veto additional jets 




Et>S GeV, 2 <| 77 |< 3.6 




A0(Jetl, Jet2) accepted range (deg) 


70 - 160 


70 - 160 


70 - 160 


A(;/)(Jet, Ifr) accepted range (deg) 


45 - 180 


45 - 180 


45 - 180 


£T(Jet2) + 9t (GeV) 


> 65 


> 85 


> 105 


ET,^.M.T (GeV) 


< 15 


< 15 


< 15 


Lepton Veto 


YES 


YES 


YES 


minimum # tracks in jet (1 1< 1) 


4 


4 


4 


JET-PROBABILITY Tagging 




1 tag ( JP < 5%) 





TABLE IL The event selection cuts for pair production of sbottom in the "low", "medium", and "high" sbottom mass regions. 
Jetl, Jet2, and Jet3 are, respectively, the first, second, and third leading jets. The cuts that are listed only under the "medium" 
column are common for all three mass ranges. 



Mass range 


Low 


Medium 


High 


(GeV/c") 


< 140 


140 - 180 


> 180 


$T (GeV) 


> 50 


> 55 


> 65 


Et Jetl (GeV) 


> 35 


> 55 


> 75 


Et Jet2 (GeV) 


> 15 


> 15 


> 35 


Et Jets (GeV) 


> 15 


> 15 


> 15 


1 77 1 : Jetl, Jet2, Jet3 




1 r?i 1< 1.2, 1 7,2 1< 1.5, 1 773 |< 2.0 




Veto additional jets 




Et>8 GeV, 2 <| 77 |< 3.6 




A(;/)(Jetl, Jet2) accepted range (deg) 


60 - 160 


50 - 160 


40 - 160 


A(f){Jet, J^t) accepted range (deg) 


45 - 180 


45 - 180 


45 - 180 


£;T(Jet2) + $T (GeV) 


> 80 


> 120 


> 160 


Eh.2M.T (GeV) 


< 15 


< 15 


< 15 


Lepton Veto 


YES 


YES 


YES 


minimum ^ tracks in jet ( t? < 1) 


4 


4 


4 


JET-PROBABILITY Tagging 




>= 1 tag ( JP < 1%) 





more. This selection reduces contributions from W^-fjets 
and Z-|-jets in which the gauge bosons decay into tau 
leptons. 

For the stop and sbottom search, the signal contains 
a large fraction of heavy-flavor jets compared to the SM 
background, which is dominated by light-flavor jets. To 
enhance the signal over SM background, we identify the 
heavy-flavor jets using the jet-probability (JP) algo- 
rithm (36j . Jets from heavy-flavor partons are charac- 
terized by secondary decays that are displaced from the 
primary vertex; thus their tracks have a large impact pa- 
rameter. Light-flavor jets appear to come from the pri- 
mary interaction and their tracks' impact parameter are 
consistent with the primary vertex (within the resolution 
of the tracking detector). The JP algorithm examines the 
impact parameter of each track from a candidate jet and 
computes a probability that the jet is a Ught/heavy jet. 
Jets from the primary (secondary) vertex are assigned a 
large (small) JP value. For the stop (sbottom) search at 



least one jet was required to have a. JP < 5% {JP < 1%). 
A looser cut is used to tag c jets as their lifetime is shorter 
compared to b jets. The efficiency to tag a fiducial c or 
6 jet increases with the transverse energy of the jet, and 
plateaus at Et ~ 80 GeV. The average efficiency to tag a 
fiducial c jet {b jet) with JP value of JP < 5% {JP < 1%) 
is-17%(- 40%) 0. 

A. Light-Flavor Background 

A significant source of background is that due to mis- 
identification of light-flavor jets as heavy flavor. In this 
case one or more light-flavor jets are tagged as a c jet 
or b jet by the JP algorithm (mis-tag). A detailed de- 
scription of the mis-tag rate is given in [36i |. Here we 
provide a brief summary of its measurement and how it 
is applied in the stop and sbottom searches. The mis-tag 
rate is measured in inclusive jet samples. The rate of 
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mis-identifying a light-flavor jet as a heavy-flavor jet is 
about ~ 1% 5%) for tagging at JP < 1% {JP < 5%). 
Uncertainties on the mis-tag rates are due to uncertain- 
ties in the contribution of long-hved A's and K's in the 
hght-flavor jets and uncertainties in the effect of particle 
interactions in the detector materials; these relative un- 
certainties are 9% for JP < 1%, and 13% for JP < 5%. 
An additional relative uncertainty on the mis-tag rate of 
6.7% (4.7%) for JP cut of JP < 1% {JP < 5%) is esti- 
mated by comparing the observed and predicted tag rates 
in different data samples [inclusive jet samples taken with 
different jet Et thresholds, and high jet multiplicity (> 4 
jets) sample]. To estimate the contribution from light- 
flavor background in the final data sample, we apply the 
measured mis-tag rate to the data sample, i.e. we mul- 
tiply each jet by the mis-tag probability, after all selec- 
tion cuts are applied except the heavy-flavor jet tagging. 
In the estimate of the heavy-flavor background contribu- 
tions with MC samples, we check for the existence of a & 
or a c parton or hadron within a cone of AR = 0.4 around 
the jet before we tag the it, to avoid double counting of 
background events due to light-flavor jets. Since a looser 
tagging requirement is used to tag the c jets in the stop 
search than to tag the b jets in the sbottom search, we 
expect the stop analysis to have a larger fraction of mis- 
tag background compared to the sbottom analysis. This 
can be seen in Tables |lll] and |lVl 



IV. DETECTION EFFICIENCY AND 
SYSTEMATIC UNCERTAINTIES 

The total detection efficiencies for the stop and sbot- 
tom signals are estimated using the pythia event genera- 
tor and the CDF detector simulation program. The sam- 
ples were generated using the CTEQSl PDFs, with the 
renormalization and factorization scale set to the mass 
of the squark in the search [SOj]. The total stop (sbot- 
tom) signal selection efficiency in the accessible mass re- 
gion varies from 0.1% to 3.4% (0.17% to 8.5%). The 
efficiency increases for higher stop (sbottom) mass and 
larger mass difference between ti (bi) and Xi- 

We have estimated the main contributions to the sys- 
tematic uncertainty on the signal acceptance and the 
SM background estimation. The uncertainty due to the 
jet energy scale for the SM background estimation is of 
the order of 10%, whereas it varies from 4% to 20% for 
the signal efficiency, depending on the stop and sbot- 
tom masses (larger uncertainty for smaller squark mass). 
The systematic uncertainty from NLO cross sections in 
the SM background varies between 8% to 13%. The sys- 
tematic uncertainty on the efficiency for tagging a c jet 
(6 jet) is 12% (8.6%) [36]. The uncertainty on the sig- 
nal acceptance due to modeling of gluon radiation from 
the initial-state or final-state partons is 5%. This is eval- 
uated from signal MC samples generated with different 
levels of initial and final state radiation. The uncertainty 
on the signal efficiency due to the PDF choice is deter- 



mined to be 2%, using the CTEQGm uncertainty PDF set. 
The uncertainty from MC statistics reaches in the most 
selective search region 50% for the SM background, and 
10% for signal. The uncertainty on the trigger efficiency 
is 5%, and the uncertainty on the luminosity of the data 
sample is 6% (s^j. The quoted uncertainties are relative 
to the estimated signal and backgrounds. 



V. RESULTS AND CONCLUSIONS 

The SM contributions and the total number of ob- 
served data events are shown in Tables IIIII and IIVI for 
the stop and sbottom searches respectively. We find 
that after applying all the selection cuts, the number of 
observed events are consistent with the number of ex- 
pected SM background events for both stop and sbottom 
searches. In these two searches, and for all the mass 
ranges, the largest source of background events is due to 
mis-identification of light-flavor jets as c or 6 jets. The 
HF multi-jet background is the second largest source of 
SM background in the "low" mass range. However, its 
contribution is largely suppressed by the tighter cuts em- 
ployed in the "medium" and "high" mass ranges. 

We studied several kinematic distributions. As an ex- 
ample Fig. [3] shows the observed I^t distributions from 
the data and the predicted I^t distributions from the SM 
background for the "high" mass stop and "medium" mass 
sbottom searches after all selection criteria are applied. 
They are in good agreement with the distributions ob- 
served in the data. No evidence for stop and sbottom 
production is observed in any of the kinematic regions 
that we have studied. 

For each of the three stop/sbottom mass ranges the 
number of observed events in the data is consistent with 
the SM expectation. An upper limit on the possible num- 
ber of signal events at 95 % C.L. using a modified fre- 
quentist approach [38*] is calculated for each mass range. 
We compare this upper limit to the prediction from the 
NLO calculation of prospino using the CTEQ6m PDFs. 
The uncertainties on the theoretical cross section arise 
from the choice of the renormalization and factorization 
scale and of PDFs. The changes induced by modifying 
by a factor of two higher or lower than its nominal 
value {Qj-f = m^^ ) , and the variations observed with 
the CTEQ6m uncertainty PDF set, result in a change of 
^ ±20% in the theoretical cross section, when combined 
in quadrature. To extract the limits on the stop and 
sbottom masses, we conservatively choose the one sigma 
lower bound of the NLO cross section associated with 
these uncertainties. 

The interpretation of the null result in the stop search 
is presented as a 95% C.L. exclusion region in the mass 
plane of m^o vs mj:_^, as shown in Fig. 31 The observed 
exclusion is smaller than the expected exclusion because 
we observe slightly more events in the data compared to 
the number of events expected from the SM processes, 
as shown in Table [TTTl The maximum exclusion value of 
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TABLE III: The number of observed data events, and the number of expected events from standard model sources in the stop 
signal region. The first uncertainty is from limited simulation statistics and the second is from systematic uncertainties. 



Mass range 


Low 


Medium 


High 


(GeV/c2) 


< 100 


100 - 120 


> 120 


Process 




Events expected 




W + jets 


11.5 ±2.4 ±2.6 


9.3 ±2.3 ±2.1 


4.0 ± 1.5 ±0.9 


Z + jets 


9.9 ±0.5 ±2.0 


7.3 ± 0.4 ± 1.5 


4.1 ±0.3 ±0.8 


Di- boson 


2.5 ±0.1 ±0.5 


2.0 ±0.1 ±0.4 


0.9 ±0.1 ±0.2 


Top 


5.2 ±0.2 ±0.8 


4.9 ±0.2 ±0.8 


3.9 ±0.2 ±0.6 


HF Multi-jet 


32.5 ±5.2 ±8.1 


17.9 ±4.0 ±4.5 


2.6 ± 1.5 ±0.6 


Mistag 


75.4 ±2.2 ±10.7 


53.5 ±2.0 ±7.6 


27.2 ±1.5 ±3.9 


Total Expected 


137 ± 6.2 ± 14.6 


94.9 ±5.0 ±9.9 


42.7 ±2.6 ±4.6 


Data 


151 


108 


43 



TABLE IV: The number of observed data events, and the number of expected events from standard model sources in the sbottom 
signal region. The first uncertainty is from limited simulation statistics and the second is from systematic uncertainties. The 
contribution from di-boson background is found to be negligible in the sbottom search. 



Mass range 


Low 


Medium 


High 


(GeV/c^) 


< 140 


140 - 180 


> 180 


Process 




Events expected 




W + jets 


5.5 ± 1.2 ± 1.2 


1.5 ±0.6 ±0.3 


0.5 ±0.4 ±0.1 


Z + jets 


6.0 ± 0.4 ± 1.1 


2.7 ±0.3 ±0.5 


1.0 ±0.2 ±0.2 


Top 


4.2 ±0.2 ±0.6 


2.9 ±0.1 ±0.4 


1.2 ±0.1 ±0.2 


HF Multi-jet 


18.8 ±4.0 ±4.7 


2.6 ± 1.5 ±0.7 


n+2.1 


Mistag 


20.5 ±0.6 ±2.3 


8.1 ±0.5 ±0.9 


2.0 ±0.3 ±0.2 


Total Expected 


55.0 ±4.2 ±5.9 


17.8 ± 1.7 ± 1.6 


4.7l?:5±0.5 


Data 


60 


18 


3 



TOjj^ is 132 GeV/c^ for m-^o = 48 GeV/c^, corresponding 
to a cross section times branching ratio limit of 2.8 pb. 
The maximum m^o excluded is 57 GeV/c^ at m^^ = 120 
GeV/c^. The reach in m^^ is limited by the integrated 
luminosity, whereas the gap between the kinematic limit 
m^j^ = TTic ± rUj^o and the excluded region is mostly due 
to the I^T requirement in the event selection. 

For the sbottom search, the interpretation of the null 
result is presented as a 95% C.L. exclusion region in the 
mass plane of m-j^o vs jti^^, as shown in Fig. \5\ The 
observed limit is larger than the expected limit at large 
sbottom mass because we observe fewer events than ex- 
pected from the SM processes in the "high" sbottom mass 
range search. The mg^ is excluded up to 193 GeV/c^ for 
m^o = 40 GeV/c^. This corresponds to a cross section 
times branching ratio limit of 0.25 pb. The exclusion fea- 
tures of this plot are similar to those for the stop search. 

In conclusion, we have performed searches for stop 
and sbottom production in proton-antiproton collisions 
at = 1.96 TeV, using 295 pb~^ of data recorded by the 
CDF experiment. In this search we assume that the stop 
(sbottom) decays exclusively into a c (b) quark and the 
lightest neutralino. The number of events that pass our 
selection for both stop and sbottom searches is consistent 
with the standard model expectation. Our 95% C.L. ex- 



clusion regions in the supersymmetric quark-ncutralino 
mass plane extend be yon d the LEP and the Tevtron's 
Run I reaches [l^, [T]| . The exclusion reach in this 
stop search is comparable with D0's latest Run II stop 
results fl6j . 
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FIG. 3: The Ifx distribution in the stop (top) and sbottom 
(bottom) signal regions for data (sohd points) compared to 
the SM background (fiUed histograms). Also shown (open 
histograms) are the expected distributions arising from stop 
pair production and decay at rrij^ — 120 GeV/c^ and m^o — 

50 GeV/c^, and sbottom pair production and decay at m^^ — 
160 GeV/c^ and m^o = 80 GeV/c^ 
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FIG. 4: The 95% C.L. exclusion region in the mass plane of 
m^o vs mi_^, assuming that stop decays exclusively into cxi'. 
For the region excluded by the present search (area inside the 
solid curve), the one sigma lower bound of the predicted NLO 
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